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OPserver: interactive online-computations of opacities and 
radiative accelerations 
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ABSTRACT 

Codes to compute mean opacities and radiative accelerations for arbitrary chemical 
mixtures using the Opacity Project recently revised data have been restructured in a 
client-server architecture and transcribed as a subroutine library. This implementation 
increases efficiency in stellar modelling where element stratification due to diffusion 
processes is depth dependent, and thus requires repeated fast opacity reestimates. 
Three user modes are provided to fit different computing environments, namely a web 
browser, a local workstation and a distributed grid. 

Key words: atomic processes - radiative transfer - stars: interior. 



1 INTRODUCTION 



Astrophysical opaci ties from the Opacity Project (OP) have 
been updated by iBadnell et al.l (|2005h to include inner- 
shell contributions and an improved frequency mesh. The 
complete data set of monochromatic opacities and a suite 
of codes to compute mean opacities and radiative accel- 
erations (OPC D.2.£l) have also been publicly released by 
ISeatonl (|2005l ) to make in-house calculations for arbitrary 
mixtures more versatile and expedient. Regarding data ac- 
curacy, there is excellent overa ll agreement between the 
OPAL dlglesias fc Rogers! 1 199ft) an d OP results a s dis - 
cuss ed by ISeaton fc Badnelll (120041) . IBadnell et all (120051 ) 
and lDelahave fc Pinsonneaultl 1 2005t ). 

Rosseland mean opacities are sensitive to both the 
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basic atomic data used and the assumed abundances of 
the chemical elements. What had been a good agree- 
ment between theory and the helioseismological data was 
found to be le s s goo d using revised solar abundances from 
lAsplund et al.l (|2005l ). The revised OP opaci ties have been 
instrumental in discussions of th at problem |Antia fc Basul 
| 2005l; iBahcall et all l2005allbll3 ; iBahcall fc Serenellil 120051 
iDelahave fc Pinsonneault 20061 ) . 

The modelling of stellar interiors, on the other hand, is 
being renewed with the solar experience. Present ( WIRE, 
MOST, CoRoT) and future (Kepler) space probes and 
the well established solar methods are giving the field 
of asteroseismology re markable growth and the guaran- 
tee of invaluable data (iMetcalfe et al.l 12004 iKurtd 120051 : 
IChristensen- Dalsgaardll2006l ). In future work on stellar mod- 
els it may be desirable to take account of revisions in abun- 
dances similar to those performed for the Sun. 

For some types of stars, models must take into ac- 
count microscopic diffusion processes, e.g. radiative levita- 
tion, gravitational settling and thermal diffusion, as they can 
affect the internal and thermal structures, the depth of the 
convection zone, p ulsations and give rise to surface abun- 
dance anomalies (jSeatonl 1 19991 ; IDelahave fc Pinsonneaultl 
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120051; iBourge fc Alecianl |2006h . As reviewed by iMichaudl 
( 20041 ). such processes are relevant in the description of 
chemically peculiar stars, horizontal-branch stars, white 
dwarfs and neutron stars, and in globular cluster age de- 
terminations from Population II turnoff stars. Furthermore, 
in order to solve the outstanding discrepancy of the at- 
mospheric Li abundance in old stars with that predicted 
in big-bang nucleosynthesis, iRichard et all (|2005Tj have pro- 
posed Li sinking deep into the star due to diffusion. This 
hy pothesis has been recently confirmed in the observations 
bv lKorn et ail (|2006l ). 

The OPCD_2.1 release includes data and codes to com- 
pute the radiative accelerations required for studies of dif- 
fusion processes. It should be noted that the radiative ac- 
celerations are summed over ionization stages and that data 
for the calculation of diffusion coefficients are calculated as- 
suming that the distribution over ionization stages of the 
diffusing ions is the same as that in the ambi ent plasma. The 
validit y of this approximation is discussed bv lGonzalez et al.l 
(I 19951) . 

In some cases, particularly when element stratification 
depends on stellar depth, calculations of mean opacities and 
radiative accelerations must be repeated at each depth point 
of the model and at each time step of the evolution, and 
thus the use of codes more efficient than those in OPCDJ2.1 
may be necessary. This becomes critical in the new dis- 
tributed computing grid environments where the network 
transfer of large volumes of data is a key issue. In the present 
work we have looked into these problems, and, as a solu- 
tion, report on the implementation of a general purpose, 
interactive server for astrophysical opacities referred to as 
OPserver. It has been installed at the Ohio Supercomputer 
Center, Columbus, Ohio, USA, from where it can be ac- 
cessed through a web pag43 or a linkable subroutine library. 
It can also be downloaded locally to be run on a stand-alone 
basis but it will demand greater computational facilities. In 
Section 2 we discuss the computational strategy of the codes 
in OPCD_2.1 followed by a description of OPserver in Sec- 
tion 3. In Section 4 we include some tests as an indication 
of its performance with a final summary in Section 5. 
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where is the mean atomic weight. The cr(u) is a weighted 
sum of the monochromatic opacity cross sections for each of 
the chemical constituents 



and is conveniently tabulated on the w-mesh 
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and Vao = v(u — » 00). The rationale behind the u-mesh is 
that it enhances fre quency resolution where F(u) is large 
|Badnell et al.ll2005l ). 



2.2 Radiative accelerations 

Similarly, the RA for a selected k element can be expressed 
as 

fJtKR'Jk 
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where fi k is its atomic weight and c the speed of light. The 
function T is given in terms of the effective temperature T e ff 
and fractional depth r/R+ of the star by 

T = itB(T cS )(R i ,/r) 2 (7) 

with 

2(nk B T) 4 



B{T) 15c^3 ' 
The dimensionless parameter 
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depends on the cross section for momentum transfer to the 
k element 



2 OPCD CODES 

We highlight here some of the key features of the codes in 
OPCD_2.1. For a chemical mixture specified by abundance 
fractions f k , they essentially compute two types of data: 
Rosseland mean opacities (RMO) and radiative accelera- 
tions (RA). 



(7^ ta = CT fc (u)[l - exp(-it)] - a k (u) 



(10) 



where a k (u) is a correction to remove the contributions of 
electron scattering and momentum transfer to the electrons. 
Both <Jk(u) and ak{u), which are hereafter referred to as the 
mono data set (~1 GB), are tabulated in equally spaced v 
intervals to facilitate accurate interpolation schemes. 



2.1 Rosseland mean opacities 

For the frequency variable 

u = hv/k B T (1) 

where k B is the Boltzmann constant, RMO are given by 
the harmonic mean of the opacity cross section a(u) of the 
mixture 



2.3 Computational strategy 

The computational strategy adopted in the OPCD_2.1 re- 
lease is depicted in the flowcharts of Figure 1 where it may 
be seen that calculations of RMO and RA are carried in 
two stages. In a time consuming Stage 1, RMO and RA 
are computed with the mixv and accv codes, respectively, 
on a representative tabulation of the complete temperature- 
electron-density (T, A^e) plane. In mixv the chemical mixture 
is specified in the input file mixv . in as 
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{X,Z,N,Z k J k } 



(11) 
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http: / / vizier . u-strasbg.fr /topbase / opserver /figl .eps 

Figure 1. Flowcharts for the computations of Rosseland mean opacities (RMO) and radiative accelerations (RA) with the codes in the 
OPCD_2.1 release. They are carried out in two stages: in a time consuming Stage 1, data are computed for the whole (T,N e ) plane 
followed by fast bicubic interpolations in Stage 2. The intermediate files mixv.xx and acc.xx enable communication between these two 
steps. 



where X and Z are the hydrogen and metal mass-fractions, 
N the number of elements, and Zk and ft are the metal nu- 
clear charges and fractional abundances. In accv, the input 
data (accv . in) are 

{N,Z k ,f h ,Zi,N x ,xs} (12) 

where now k runs over the N elements of the mixture, and 
Zi and Xj are respectively the nuclear charge and N x abun- 
dance multipliers of the test i element. Input data formats 
in either mixv.in or accv. in give the user flexible control 
over chemical mixture specifications. 

As shown in Figure 1, the intermediate output files 
mixv.xx (~85 KB) containing 

{p,Kn}(T,N e ) , (13) 
where p is the mass-density, and acc.xx (~470 KB) with 

{«"^.£}<™.»> <»' 

are written to disk. They are then respectively read by 
the codes opfit and accfit in Stage 2 for fast bicubic 
interpolations of RMO and RA on stellar depth profiles 
{T, p, specified by the user in the opfit. in and 

accfit. in input files. The final output files are opfit.xx 
containing 

/ dlogKR 91ogK R l 

SlogKR,^- — W (15) 

[ ologT a log p J 

and accfit. xx with 

{log«: R ,log7,logg ra d}(i,Xj) • (16) 

In this computational approach, performance is mainly 
limited by the summation in equation @ which implies disk 
reading the mono data set; for instance, in mixv it takes 
up to ~90% of the total elapsed time. OPCD_2.1 also in- 
cludes other codes such as mx and ax which respectively 
compute RMO and RA for a star depth profile. The chemi- 
cal mixture can be fully varied at each depth point using the 
specifications in equations (|llH12p . the RMO and RA being 
obtained in a one-step process using bicubic interpolations 
without splines. These methods are thus suita ble for cases 
with multi-mixture depth profiles l|Seatonll2005l ) . Further de- 
tails of all the OPCD codes are contained in the reference 
manuajf]. 



3 OPSERVER 

In OPserver the computational capabilities of the codes in 
OPCD_2.1 are greatly enhanced by the following innovative 
adaptations. 

3 http : //opacities . osc . edu/publi/OPCD .pdf 



(i) The codes are restructured within a client -server net- 
work architecture whereby the time consuming Stage 1 is 
performed on a powerful processor while the fast Stage 2 is 
moved to the client, e.g. a web server or a user workstation. 
In this arrangement performance could be affected by the 
client-server transfer of the mixv.xx and acc.xx intermedi- 
ate files, but since they are never larger than 0.5 MB, it is 
not expected to be a deterrent with present-day bandwidths. 
In a local installation where both the client and server reside 
on the same machine, communication is managed through 
shared buffers in main memory; i.e. the corresponding data 
in mixv.xx and acc.xx are not written to disk. 

(ii) The codes are transcribed as a subroutine library — 
to be referred to hereafter as the OPlibrary — which can be 
linked by the user stellar modelling code for recurrent sub- 
routine calls that avoid data writing on disk. That is, the in- 
put data in the mixv. in, accv. in, opfit . in and accfit . in 
files and the output tables in the opfit.xx and accfit.xx 
files (see Figure 1) are now handled as subroutine param- 
eters while the intermediate mixv.xx and acc.xx files are 
passed via shared main-memory buffers. Chemical mixtures 
are again specified with the formats of equations (|llH12fl 
which allow full variation at each depth point in a single 
subroutine call. 

(iii) RMO/RA are computed with the complete mono data 
set always loaded in main memory thus avoiding lengthy and 
repeated disk readings. This is achieved by implementing 
OPserver on a dedicated server where mono is permanently 
resident in RAM, or in the case of a local installation, by 
disk-reading it once at the outset of a modelling calculation. 

(iv) When accessing the remote server, client data re- 
quests are addressed through the HTTP protocol, i.e. in 
terms of a Uniform Resource Locator (URL). This allows 
data fetching from the central facility through an interactive 
web page or a network access subroutine, the latter being 
particularly suitable for a stellar model code that is to be 
run in a distributed grid environment. 

(v) The do-loop that computes the summation of equa- 
tion ([3]) has been parallelized in OpenMP which provides 
a simple, scalable and portable scheme for shared-memory 
platforms. 

As shown in Figure 2, the current OPserver enterprise is 
implemented as a client -server model at the Ohio Supercom- 
puter Center (OSC). The web server communicates with the 
supercomputer via a socket interface. Earlier versions were 
developed on an SGI Origin2000 server with the PowerFor- 
tran parallelizing compiler. The current version runs on a 
Linux system with Fortran OpenMP directives. OPserver 
offers three user modes with full functionality except when 
otherwise indicated in the following description. 

Mode A In this mode OPserver is set up locally on a 
stand-alone basis (see Figure 2). The facilities of the OSC are 
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Figure 2. OPserver enterprise showing the web-server-supercomputer tandem at the Ohio Supercomputer Center (OSC) and the three 
available user modes. (A) The OPlibrary and monochromatic opacities (mono) are downloaded locally and linked to the user modelling 
code such that RMO/RA are computed locally. (B) The OPlibrary is downloaded locally and linked to the modelling code but RMO/RA 
are computed remotely at the OSC. (C) RMO/RA computations at the OSC through a web client. 



not used. A new OPCD release (OPCD_3.£|) is downloaded, 
followed by (i) installation of both the OPlibrary and the 
mono data set and (ii) linking of the OPlibrary to the user 
modelling code. Computations of RMO/RA are preceded 
by the reading of the complete mono data set from disk and 
therefore requires at least 1 GB of RAM. 

Mode B In this mode, the OPlibrary is downloaded, in- 
stalled and linked to the user code, but Stage 1 is performed 
remotely at the OSC (see Figure 2). This option has been 
customized for stellar modelling in a distributed grid envi- 
ronment that would otherwise imply (i.e. Mode A) the net- 
work transfer, installation and disk-reading of the mono data 
set at runtime. It is also practical when local computer capa- 
bilities (RAM and/or disk space) are limited. The functions 
provided by the mx and ax codes have not been implemented. 

Mode C In this mode RMO/RA computations at the 
OSC are requested through an interactive web pag^fl which 
allows both Stage 1 and Stage 2 to be carried out re- 
motely or, alternatively, Stage 2 locally by downloading the 
mixv.xx/acc.xx intermediate files (see Figure. 1) with the 
browser for further processing with local opfit/accfit ex- 
ecutables. 



4 TESTS 

OPserver benchmarks were initially carried out on an SGI 
Origin2000 multiprocessor at the OSC with an earlier re- 
lease of OPCD. For the standard S92 mixture ijSeaton et al.l 
1994), the mixv code took up to 140 s to compute the 
mixv.xx file, of which 126 s were dedicated to disk-reading 
and 14 s to the actual computing of the mean opacities. 
OPserver took on average 12.0 ± 0.5 s to compute mixv.xx 
which was not written to disk unless requested. In Fig- 
ure 3 we show the acceleration obtained on the Origin2000 
through parallelization where the calculation of mean opac- 
ities is reduced to 2 s with 8 processors. Further significant 
acceleration is prevented by data transfer overheads. 

On more recent workstations, the local performances of 
the codes in OPCD_2.1 and OPserver depend on processor 
speed and RAM and cache sizes. For instance, on a Pow- 
erMac G5 (PowerPC 970fx processor at 2.0 GHz, 1GB of 
RAM and L2 cache of 512 KB) the first time mixv is run it 
takes for a single S92 mixture 103.8 s to compute the RMO, 
but on subsequent runs the elapsed time is reduced to an av- 
erage of 28.2 ±0.2 s. Similarly, OPserver takes 103.3 s which 
is then reduced to 31.4 ± 0.4 s on subsequent runs. Once 
the mono data set is loaded in RAM by OPserver (Mode 
A), calculations of RMO for a single S92 mixture only take 
5.29 ± 0.02 s and 6.13 ± 0.01 s for RA for the test element 
Ar. In Mode B, where Stage 1 is carried out remotely at the 
OSC and the mixv.xx and accv.xx files are transferred at 

4 http : //cdsweb . u-strasbg . f r/topbase/op . html 
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the relatively low rate of 1.88 KB/s, computations of RMO 
and RA take 5.9±0.1 s and 9.3±0.3 s, respectively. The no- 
ticeable longer time taken for the latter is due to the transfer 
time taken for the larger accv.xx file. 



5 SUMMARY 

Rosseland mean opacities and radiative accelerations can be 
computed from OP data in any one of the following ways. 

(i) Download th e original OPCD_2.1 package as described 
bv lSeatonl (|2005l ) and perform all calculations locally. 

(ii) Mode A, download the upgraded OPCD_3.3 package, in- 
stall OPserver and perform all calculations locally by link- 
ing the subroutines in the OPlibrary. Calculations with 
OPserver are more efficient but require large local computer 
memory. 

(iii) Mode B, as Mode A but with Stage 1 performed re- 
motely at the OSC. Mode B is convenient if fast calculations 
are required but local computer memory is limited or when 
stellar modelling is to be carried out in a grid environment. 

(iv) Mode C, perform all calculations remotely at the OSC 
through an interactive web page whereby files are down- 
loaded locally with the browser. 
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